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The glitch activity of young pulsars arises from the exchange of angular momentum between the 
crust and the interior of the star. Recently, it was inferred that the moment of inertia of the crust of 
a neutron star is not sufficient to explain the observed glitches. Such estimates are presumed in the 
Einstein's general relativity in describing the hydrostatic equilibrium of neutron stars. The crust of 
the neutron star has a space-time curvature of 14 orders of magnitude larger than that probed in solar 
system tests. This makes gravity the weakest constrained physics input in the crust related processes. 

We calculate the ratio of the crustal to the total moment of inertia of neutron stars in the scalar-tensor 
theory of gravity and the non-perturbative f{TV} = TZ + alZ^ gravity. We find for the former that 
the crust to core ratio of the moment of inertia does not change significantly from what is inferred 
in general relativity. For the latter we find that the ratio increases significantly from what is inferred 
in general relativity in the case of high mass objects. Our results suggest that the glitch activity of 
pulsars may be used to probe gravity models, although the gravity models explored in this work are 
not appropriate candidates. 

PACS numbers: 


I. INTRODUCTION 


Rotationally powered pulsars are strongly magnetized rapidly rotating neutron stars spinning-down under the 
magnetic dipole radiation torques More than 2000 pulsar have been discovered to date [see0, for a review]. 

Some of these objects, predominantly the younger ones, show occasional and abrupt spin-up events called glitches 
Si of the size AD./Cl ~ 10 ^ — 10 ^ where D is the rotational angular velocity of the neutron star [see i E3- 
for observational data]. Of fhe more than 2000 pulsars about 160 were observed to glitch with a total of ~ 400 
glitches. The glitches are generally accompanied by a sudden change of the spin-down rate | AQ/ Oj lO-'^- 10-2. 
These glitches tend to relax back to their original pre-glitch state values. The long time-scale (days to years) of the 
post-glitch relaxation implies the presence of a superfluid interior iH [seeEi for a review]. 

The small glifches, as observed in the Crab pulsar, can be explained with the starquake model lll3ll in which the 
solid crust resisting the reduction of the oblateness of the figure of the spun-down star breaks at a critical strain. 
The starquake model predicts a time interval between the large glitches much longer than observed e.g. in the Vela 
pulsar [141] . More importantly the energy dissipation in such a large glitch would lead to an observable change in the 
surface temperature exceeding the bounds IlSll . Hence the starquake model can not be the mechanism underlying 
the Vela-like large glitches which comprise the majority of the observed glitches. 

Such large glitches are explained by the exchange of angular momentum between the superfluid component per¬ 
meating the inner crust of a neufron sfar and fhe relatively slower rotating outer crust. The spin-down torques due to 
the magnetic dipole radiation act on the outer crust and slow it down. The neutron superfluid permeating the inner 
crust follows the spin-down of the outer crust as they are weakly coupled by mutual friction forces ifl^ . The angular 


momentum of the superfluid is carried by an array of quantized vortices ml Il8i] and its slow-down is achieved by 


the radial motion of the vortices away from fhe rotation axis Ill9il . The pinning of the vortices to the nuclei precludes 
their radial motion preventing the spin-down of the neutron superfluid, and results in a lag in which the superfluid 
rotates more rapidly than the crust lEilli- This lag induces a Magnus force acting on the vortices and produces 
a stress on the crust. When the lag reaches a critical threshold the vortices are suddenly unpinned in an avalanche 
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transferring angular momentum from the rapidly rotating neutron superfluid to the crust causing its observed spin- 
up 13- 3]- The post-glitch relaxation may result from repinning of the unpinned vortices to new pinning sites and 
restart to "creep" in steady state 11911 . 

Superfluid neutrons of the inner crust, being in Bloch states of the crust lattice, have an effective mass larger than 


the bare neutron mass 11211,12211 . This so called 'entrainment' results with a less mobile crustal superfluid because of 


the reduced ratio of the effectively free neutron superfluid that can store and exchange angular momentum with the 
lattice |23-31- Iri order to address this problem authors of ref. f3l invoked the contribution of the toroidal magnetic 
flux region in the outer core of the star to the moment of inertia of the crust. Piekarewicz et al. 12711 , on the other 
hand, considered nuclear physics processes that could lead to a thicker crust. 

The arguments leading to the problem of the deficiency of the crustal moment of inertia for the glitch phenomena 
and explanations for addressing the problem rely on the calculations of the structure of neutron stars within the 
framework of Einstein's general relativity (GR). Although GR is tested to a precision of 10“^ in the solar system tests 
and found to be in excellent agreement with measurements [see [3 for a review], the presumption of the validity of 
GR for describing the hydrostatic equilibrium of neutron stars is a remarkable extrapolation, 5 orders of magnitude 
in compactness GM/c^R) and 14 orders of magnitude in curvature GM/c^R^) 13]. The equation of state (EoS) 
of the outer crust llsol] and the inner cru st llsill are fairly well understood [seejs^- for a review] where the gravity of the 
neutron star attains its strongest value 13311 . Gravity is then the least constrained physics input for the glitch activity 
in the crust and, thus, such phenomena may be used to constrain gravity models. The purpose of this paper is to 
draw attention to the role of gravity in the crust by calculating, for the first time, the fraction of the crustal moment 
of inertia in different models of gravity. 


II. METHOD 

n this section we briefly present the basic equations for f{TZ) and scalar-tensor theory (STT) of gravity, which are 
presented more thoroughly in 

The mathematical equivalence between the f{TZ) theories and STT is exploited in this work. Eor mathematical 
simplicity it is useful to study f{TZ) and scalar-tensor theories in the so-called Einstein frame, however all final 
results we present here are transformed to the Jordan frame (the physical frame). The explicit form of the function 
f{TZ) we adopt is f{TZ) = TZ + oR}. 

The filZ) theory action in Jordan frame is: 

S = J \/ ~S T ) T Smarter {S}tv/X) - (1) 

where Smarter is the action of the matter fields collectively denoted by x- The general form of the STT action in Jordan 
frame is the following: 

s = J (F(cl>)7^ - ( 2 ) 

2L1(<I’)) -|- Smarter {SfiV/X') - 

where O is the gravitational scalar field. F(‘I>),Z(<I>), and U{<i>) are functions that govern the dynamics of the 
gravitational scalar. In order the gravitation to carry positive energy F(T>) > 0, and the nonnegativity of the scalar 
field implies 2F(<I>)Z(<I>) -|- 3{dF (<3>) /> 0. We are using the mathematical equivalence between f{TZ) theories 
and STT theories to rewrite the action (|T]| like an STT action. For mathematical simplicity we rewrite both actions in 
Einstein frame, using a conformal transformation = F{^)giiv and the new scalar field q> defined by 

(dcpY^ _3 ( dlnFj^W I Z(0) 

\d^) 4 ^ do ) 2F{^)' ^ ’ 


We also define 
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A((p) = F-2(0), 2V{(p) = !J(<1>)F-2(0). (4) 

Now we have one Einstein frame STT action describing both theories: 

^ = 1^ / (5) 

-V(cp)) + 

^matter 

In (|5} all quantities marked with asterisk are with respect to the Einstein frame metric g*^^. 

We consider stationary and axisymmetric space-time as well as stationary and axisymmetric fluid and scalar field 
configuration. The metric describing this space-time in Einstein frame has the explicit form 

dsl = —+ r^{d6^ -|- sin^ 9d^^) — 

2co{r, 6)sin^ 6dMt, (6) 

where we kept only first-order terms in the angular velocity H (slow rotation approximation). 

The field equations in Einstein frame for both theories are as follows 


d^cp 

dr^ 


1 d r 


r{l — e = 8nGA*{cp)p + e 


-2A 


dr J 




2 -2A 


-e 

r 




dr r 

( * “ f - + 7) * = (p-3p)^'‘ + 


dr J 




,2A I 2A 

t' / 


4 dcp 


dp 


d(p 


dcp 


dr 


e‘1’ 


-A 


r^ dr 




* = 1' 

dcd{r) 
dr 


= 16TzGA'^(cp){p -h p)ci>{r), 


where we have defined 


Oi{(p) 


dlixAicp) 

-;- and CO = a — cv 


(7) 

( 8 ) 
(9) 

( 10 ) 

( 11 ) 


( 12 ) 


Here p and p are the pressure and energy density in the Einstein frame and they are cormected to the Jordan frame 
ones via p^ = A^{(p)p and p^ = A^{(p)p, respectively. 

The above system of equations, combined with the equation of state for the star matter and appropriate boundary 
conditions, describe the interior and the exterior of the neutron star. We obtain the exterior solutions by imposing 
the conditions p = p = 0. 

The boimdary conditions that we employ are the natural ones, namely at the center of the star we have p(0) = 
Pc/ A(0) = 0, ^ (0) = 0, and at infinity we have lim^_>.co (p{r) = 0,lim,._>oo (p{r) = 0, lim^-^oo d) = O. The coordinate 
radius of the star is determined by the condition p{r^) = 0 and the physical radius is given by Rg = A[^(rg)]rg. 

Imposing a different form of the conformal factor A{cp) and the potential V(cp) we can examine f{TZ) theory and 
STT with the same set of equations. In the case of 72.-squared gravity they have the form 

AHcp)=e-T3‘P A = y(<p) = l(^l-e-^)", (13) 

and for STT with no potential we adopt 


A{q)) = a = ficp V{cp)=0. 


(14) 
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The moment of inertia f of a neutron star is defined in the standard way 



where / is the angular momentum. For convenience in the numerical calculations, the integral expression 


1 = 


8tzG rs 

Jo 


A^{cp){p + p)e^ 



dr 


(15) 


(16) 


will be used. 

From now on we shall use the dimensionless parameter a —t a/Rg and the dimensionless moment of inertia 
I —> I/MqRq where Mq is the solar mass and Rq is one half fhe solar gravifational radius Rq = 1.47664 km. 


III. NUMERICAL RESULTS 

We solved, numerically, the hydrostatic equilibrium equations to investigate the behaviour of fhe ratio of /crust 
fo fhe total moment of inertia I of neufron stars for four realistic EoS', namely FPS, SLy, APR4 and MSI. The first 
one, FPS, is a soft EoS with maximal mass less than the observed maximum mass of 2 M© SHi. SLy and APR4 
have maximal masses slightly above 2 Mq, and the last one, MSI, is the stiffest EoS—it gives highest masses and 
biggest radii. Thus a wide range of stiffness is covered by fhese choices of EoS'. In Eig. [T]we present the mass of 
radius relation for all mentioned EOS. The results are presented for GR, and for both alternative theories we are 
investigating. We define fhe crusf of fhe neufron sfar as fhe layer wifh density lower than 1.5 X 10^“^ g cm ^ 

We examined different values for this critical density and saw that the precise value does not affect our results and 
conclusions. 


A. Scalar-tensor theory of gravity with vanishing potential 


We first investigate, numerically, the ratio /crust / / for neufron stars moment of inerfia in scalar-tensor theories 
of gravity with vanishing potential js^. The relativistic stellar models in scalar-tensor gravity can significa ntly 
deviate from fhe predictions of general relafivify Spontaneous scalarization is induced for ^ < —4.35 jilH . 
Yet observational studies on pulsar timing provides the constraint ^ —4.5 |28, l^l • The value of the parameter 
we adopt is fi = —4.5, which is consistent with the observational results. We did the calculations for all EoS' and 
compared the results with the results within GR. 

In Eig. |2]we present the ratio /cmst/ / as a function of the mass of the neutron star. The GR case is represented by 
a continuous black line and the STT one is in (red) dashed line. The chosen interval of masses is consistenf with the 
observations, but for representing the results for EoS MSI we presenf a bigger mass range because the STT deviations 
from GR occur only for masses close to the maximal ones for this EoS. Eor all EoS' we obtain qualitatively the same 
result: the /crust// ratio decreases rapidly with the mass of the neutron star. Eor models with M = O.SM© the ratio is 
between 0.1 and 0.15 (not presented in the graphs), and for models with maximal mass it is rmder 0.03 for all EoS'. 
The insets shows the discrepancy between GR and STT solutions. The discrepancy we define as 


A = 


■^(STT) - ^^(GR) 


Icrm(GR) 


(17) 


The difference is much smaller than 3% which is expected because the deviations from pure GR for /3 = —4.5 for 
non-rofating (slowly rofating) neufron stars are not large, for all EoS'. Eor the case of rapidly rofating neutrons 
stars, however, deviations in the moment of inerfia increase 11421] , therefore fhe differences for /crust// also increase 
as expected. As a result the deficiency in fhe crustal moment of inertia carmot be explained by the concrete STT 
examined in this paper. 
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FIG. 1: Mass of radius relation for all examined EoS'. The case of GR is represented by a solid black line, and the STT with 
/5 = —4.5 is given by black line with symbols, and the Ti? gravity results are represented by lines with different patterns and 
different colours. 


B. 72.-squared gravity 

We continue by investigating numerically the behaviour of the crust to the total moment of inertia ratio in 72.^ 
gravity. The observational constraints for the parameter a are the following - a ^ 5 or in dimensionless 

units fl < 10®. 

We find that unlike STT the TZ^ gravity leads to serious deviations from pure general relativity as presented in Fig. 
|3]for a range of values of a. For all EoS', we observe qualitatively the same behaviour: for low (high) masses the ratio 
is the highest (lowest) for GR, the cross-over occurring at around M ~ l.SM©. 

In the insets we present the discrepancy between pure GR and the results in TZ^ gravity for different values of the 
parameter a. The discrepancy we define as 

^ ^ W(;(7^))-W(GR) 

i^(GR) ' ^ 

and the presented results are for models with large masses up to the maximal ones in GR for the corresponding EoS. 
The deviation from GR with the increase of a is clearly visible. The insets show that the increase of the crustal to the 
total moment of inertia ratio in the TZ^ gravity for the maximal GR mass compared to the pure GR case is more than 
50%. 

Eor TZ^ gravity the problem with the deficiency of the crustal moment of inertia is partially solved—for models 
with masses higher than, roughly, l.SM©, the ratio increases, with a compared to GR. Eor models with small masses, 
however, the ratio decreases with a, so the moment of inertia deficiency problem gets even worse. 
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FIG. 2: Crust to the total moment of inertia ratio as a function of mass for all EoS'. The case of GR is represented by a solid black 
line, and the STT with /5 = —4.5 - by a dashed (red) line. The inset shows the discrepancy between pure GR and the results in 
STT. 


IV. CONCLUSIONS 

Given the identified deficiency of fraction of crusfal moment of inertia to the total moment of inertia of neufron 
stars (fcmst/ 1) —a critical parameter in the models addressing the observed glitch phenomena—we have investigated 
this ratio in scalar tensor theory of gravify with the conformal facfor A(^) = and 72.-squared gravify model 
(/(72.) = TZ + aR?-) to see whether the free parameters of these theories could account for fhe discrepancy. 

In fhe case of STT we have foimd thaf fcmst/f changes much less fhan 3 % for ^ = —4.5, fhe highest value 
compatible with solar system tests. This is an expected result, if we fake into account the fact that the differences for 
all neufron sfar paramefers are relatively small. This implies that the deficiency of fcmst/f can not be accounted for 
by referring fo STT. 

In fhe case of /^.-squared gravify model we have found thaf for masses M < l.SM© the ratio fcmst/ f is even smaller 
than its value in GR with increasing values of a which renders the crustal moment of inertia deficiency problem even 
more acufe. For larger masses, however, we observe fhat fcmst/f is greater than its value in GR by some 20% to 60% 
depending on the value of fhe parameter a and the mass. This result suggests that the crustal moment of inertia 
deficiency could be accounfed only if all observed glifching pulsars have masses larger fhan M > 1 .SM©. It should, 
however, be emphasized that young glitching pulsars have measured masses 1.4 M© at which the difference in GR 
and /^.-squared gravify is fhe smallest. 

We then conclude that the two gravity models we have considered in this paper are not candidates for addressing 
fhe crusfal deficiency problem. This does nof exclude thaf fhe problem could be accounfed for in differenf models 
of gravify. Nor we can exclude an astrophysical explanation of the problem such as the one given in Ref. 
Additional studies of alfernafive theories should be done as to allow more general conclusions to be made. For 
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FIG. 3: Crust to the total moment of inertia ratio as a function of fhe mass of the neutron star in the case of 72. + aTi? gravity. GR 
solutions are in continuous black line and different colours and patterns correspond to different values of a. The insef shows fhe 
discrepancy between pure GR and the results in 72^ gravity for different values of the parameter a. 


example, the investigation of the problem in the context of STT with a massive scalar field is in progress and the 
results will be presented elsewhere. 
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